00 ' 

Envelope burning over-luminosity: a challenge to synthetic 



in 

(N 



> 
(N 



TP-AGB models 



Paola Marigo 



1,2 



m 

' ^ Max-Planck Institut fiir Astrophysik, Karl-Schwarzschild-Str. 1, D-87540 Garching bei Miinchen, Germany 
OO . ^ Department of Astronomy, University of Padova, Vicolo dell'Osservatorio 5, 35122 Padova, Italy 



submitted to Astronomy & Astrophysics 



o 

Abstract. Until recently synthetic AGB models had not taken into account the break-down of the 
core mass-luminosity (Mc — L) relation due to the occurrence of envelope burning in the most massive 
^ . {MZ 3.5Mq) and luminous (Mboi ^ -6) stars. 

H I Marigo et al. (1998) made the first attempt to consistently include the related over-luminosity effect (i.e. 

above the Mc — L relation) in synthetic TP-AGB calculations. The method couples complete envelope 
integrations with analytical prescriptions, these latter being presently updated with the highly detailed 
relations by Wagenhuber & Groenewegen (1998). 

In this paper the reliability of the solution scheme is tested by comparison with the results of com- 
plete evolutionary calculations for a TMq AGB star undergoing envelope burning (e.g. Blocker & 
Schonberner 1991). 

Indeed, the method proves to be valid as it is able to reproduce with remarkable accuracy several 
evolutionary features of the 7Mq star (e.g. rate of brightening, luminosity evolution as a function of 
the core mass and envelope mass for different mass-loss prescriptions) as predicted by full AGB models. 
Basing on the new solution method, we present extensive synthetic TP-AGB calculations for stars with 
initial masses of 3.5, 4.0, 4.5, and 5.0 Mq, and three choices of the initial metallicity, i.e. Z = 0.019, 
Z = 0.008, and Z = 0.004. Three values of the mixing-length parameter are used, i.e. a = 1.68, 2.0, 2.5. 
We investigate the dependence of envelope burning on such stellar parameters (M, Z, and a). The com- 
parison between different cases gives hints on the interplay between envelope burning over-luminosity 
and mass loss, and related effects on TP-AGB lifetimes. 
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1. Introduction 

Synthetic AGB models may provide a powerful tool of investigation (Renzini & Voli 1981; Groenewegen & de Jong 
1993; Marigo at al. 1996a). Evolutionary calculations can be easily carried out over the whole mass range of AGB 
stars and for various metallicities, and the results tested according to the adoption of various input prescriptions 
(e.g. analytical relations derived from full AGB models, mass-loss laws, parameters of the third dredge-up). In this 
way, it is possible to explore the sensitiveness of theoretical predictions to different physical assumptions and, at the 
same time, to readily get an overall picture of aspects related both to single star evolution (e.g. location on the H-R 
diagram, maximum AGB luminosities, stellar lifetimes, initial mass-final mass relation, changes in the surface chemical 
composition, stellar yields) and to integrated properties of the stellar aggregates in which AGB stars are present (e.g. 
luminosity functions of oxygen- and carbon-rich AGB stars, contribution of evolved stars to the integrated light and 
to the chemical enrichment of the host galaxy). 

Of course, the computational agility and flexibility typical of the synthetic approach are paid with a certain loss 
of details if compared to complete AGB models. However, in order to get reliable results from synthetic analyses the 
degree of approximation in treating the physical processes must be good enough so that no essential feature is missed. 

Actually, a clear point of inadequacy of synthetic AGB models has so far concerned the treatment of envelope 
burning in the most massive stars (M ^ 3.5Mq). In brief, the erroneous assumption common to AGB analyses 
performed either with the aid of convcctive envelope models (Scalo ct al. 1975; Renzini & Voli 1981; Marigo ct al. 
1996a), or purely analytical prescriptions (Groenewegen & de Jong 1993) is that the quiescent surface luminosity of 
a TP-AGB star experiencing envelope burning still obeys the standard Mc — L relation, as in the case of lower mass 
stars (M < 3.5Mq). 

On the contrary, over the years complete AGB calculations have clearly indicated that nuclear burning in the hottest 
convective envelope layers of the most massive TP-AGB stars may be a considerable energy source, making these stars 
depart significantly from the Mc — L relation towards higher luminosities (Blocker & Schonbcrncr 1991 (hereinafter 
also BS91); Lattanzio 1992; Boothroyd & Sackmann 1992; Vassiliadis & Wood 1993; Blocker 1995 (hereinafter also 
B95); Wagenhuber 1996). 

Indeed, the recent discovery of the over-luminosity produced by envelope burning above the Mc — L relation has 
notably changed and improved our current understanding of the AGB phase, setting important implications. 

In principle, the luminosity increase of a TP-AGB star is not bounded by the classical AGB luminosity limit of 
Afboi ^ ^7.1, predicted by the Paczynski (1970) Mc — L relation for a mass of the degenerate core equal to the 
Chandrasekhar critical value of ~ IAMq (BS91). For this reason, the observation of AGB stars with luminosities close 
to Mboi ~ — 7.1 does not imply that their core masses are close to IAMq, a circumstance that has been until recently 
considered as a supporting indirect evidence for the occurrence of Type I1/2 Supernovae, consequent to disruptive 
carbon ignition (see for example Wood et al. 1983). 

Moreover, the maximum quiescent luminosity, Lmax, attained by an AGB star with envelope burning does not 
coincide with the final luminosity at the so-called AGB tip, as in the case of lower mass stars following the core 
mass-luminosity relation. In higher mass stars the maximum of the luminosity occurs, in practice, at the onset of the 
superwind, and hence before the AGB tip is reached (when the envelope has been almost completely ejected and the 
core mass- luminosity relation is recovered) . 

It follows that at high luminosities, say Mboi < —6.5, the core-mass luminosity relation cannot be longer employed 
in combination with the initial mass-final mass (Mi — Mf) relation to estimate the age of a coeval system on the base 
of its brightest AGB stars (usually inferred through the steps: imax — > Mf — > Mi — > age; see for instance the review 
by Iben k Renzini (1983)). 

The over-luminosity effect is also expected to intensify the mass-loss process suffered by stars with envelope burning, 
so as to possibly anticipate the onset of the superwind regime. Then, this would result in a reduction of the TP-AGB 
lifetimes and hence of the remnant white dwarf masses. The latter consequence may concur to solve the long-standing 
problem related to the excess of white dwarfs more massive than ~ 0.7Mq, as predicted by synthetic AGB models 
(see Bragaglia et al. (1995) for a punctual discussion of this point). 

In this context, Marigo et al. (1998) first pointed out a possible solution scheme to overcome the limit of synthetic 
models. This paper aims at verifying the validity of the original treatment of envelope burning - coupling the use of 
analytical relationships with complete envelope integrations - by testing its capability of reproducing the results from 
full AGB calculations. 

The general organisation of the paper is as follows. In Sect. 2 the rationale underlying envelope integrations is 
briefly recalled, together with the main analytical prescriptions. In Sect. 3 the validity of the method is checked by 
comparison with full evolutionary calculations for a 7M0 star performed by BS91 and B95. Section 4 presents the 
results of synthetic TP-AGB evolutionary calculations for stars with initial masses in the range 3.5M0 < M < S.OMq, 
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and initial mctallicitics Z ~ 0.019, Z — 0.008, and Z — 0.004. The sensitiveness of envelope burning to stellar mass, 
metallicity, and mixing-length parameter is discussed in Sect. 5. Some relevant quantities characterising the evolution 
of the model stars are tabulated in Table 2. Finally, Sect. 6 contains some concluding remarks and illustrates the 
intents of future works. 



2. Outline of the method 

The reader is referred to the work by Marigo et al. (1998) for a detailed description of the method developed to 
calculate the energy contribution from envelope burning to the stellar luminosity. Let us herein just summarise the 
basic points. 

Given the total stellar mass M, the core mass Mc, and the chemical composition of the convective envelope at each 
time during the quiescent intcr-flash periods, the surface luminosity L is singled out by means of envelope integrations, 
provided that proper boundary conditions are fulfilled. 

In this way, the erroneous assumption L = ~ where Lm^ corresponds to the luminosity predicted by the 
Mc — L relation for a given core mass — is abandoned. Moreover, the usual prescription, fixing Lr = constant = L 
throughout the envelope, is replaced with the equation of energy balance under static approximation (i.e. the entropy 
term —TdS/dt is neglected): 

so that the complete set of the stellar structure equations must be integrated. 

In general, the determination of the imknown fmictions r, P^, T^, Lr across the envelope requires to specify four 
boundary conditions. For any given pair of envelope parameters, L and Teg, two conditions naturally derive from the 
integration of the photospheric equations for T and P down to the bottom of the photosphere (Kippenhahn et al. 
1967). Then, since both L and Tefj are actually free parameters in our static envelope model, two more boundary 
conditions must be fixed. 

To this aim, we express the quiescent surface luminosity L as: 

L = Lq + Lhb — + Lu + Leb (2) 

where Lq represents the rate of energy generation due to the gravitational contraction of the core: ifjc is the small 
energy contribution from the He-burning shell; is the rate of energy loss via neutrinos; ijj ^md Leb refer to the 
rate of energy production by hydrogen burning in radiative and convective conditions, respectively. 

We remark that the gravitational contribution is included only in the energetic budget via the term Lq, but 
according to the static approximation used in our model, we do not actually take into account the possible contraction 
and/or expansion of the structure during the evolution. 

It is worth noticing that the energy contributions indicated in the right hand-side of Eq. (2) are produced within 
distinct regions of the star. Specifically, {Lq + Lue ~ L^) represents the net rate of energy outflow from the core, 
this latter being commonly defined as the stellar interior below the H-He discontinuity. Beyond the core, energy is 
produced by nuclear burning of hydrogen, at a total rate given by the the sum of the two terms (Ljj + Leb)- It results 
that Leb = in low-mass TP-AGB stars (M ^ 3.5Mq) complying with the Mc — L relation (i.e. without envelope 
burning). Differently, the term Leb can significantly contribute to the energy budget of more massive TP-AGB stars 
{M > 3.5Mq), as the base of the convective envelope penetrates into the H-burning shell. 

In virtue of the site separation of the energy sources, it is convenient to set the two boundary conditions in 
question by specifying the local values of the luminosity Lr at two suitable transition points, provided that Eq. (2) is 
satisfied. Denoting by i?corc and i?conv the radial coordinates of the bottom of the H-burning shell (i.e. where hydrogen 
abundance is zero, below which the mass Mc is contained) and the base of the convective envelope, respectively, we 
can write: 

L{r = Rcore) = Lq + Liie - (3) 

L{r = Rconv) = L{r = Rcore) + (4) 

with the right hand-side members of both equations being known functions of the core mass, envelope mass, and 
chemical composition (see Sect. 2.1 for the adopted prescriptions). 

Then, Eqs. (3) and (4) together with the photospheric conditions provide the four boundary constraints necessary 
to determine the entire envelope structure. Numerical integrations of the envelope are performed with a very fine 
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mass resolution, the width of the innermost shells (where the structural gradients become extremely steep) typically 
amounting to 10"'^ - 1Q~^Mq. 

The solution model yields, in particular, the quantity Leb, i-e. the energy produced within the convective envelope, 
and the {L, Tes) pair, i.e. the current location on the H-R diagram, without invoking further external assumptions. 

2.1. Prescriptions for (Lq + Lue — L^) and Lu 

The boundary conditions expressed by Eqs. (3) and (4) imply the knowledge of the luminosity contributions: 

— (Lg + — Li,) from the core 

— Lfj from radiative hydrogen burning 

To this aim, we adopt the analytical formulae describing the light curves of TP-AGB stars, presented by Wagen- 
huber & Groenewegen (1998). These prescriptions are a high accuracy reproduction of the results from extensive grids 
of complete evolutionary calculations carried out by Wagcnhubcr (1996) for stars with initial masses in the range 
0.8Mq < M < 7.0Mq, and metallicities Z = 0.0001, Z = 0.008, and Z = 0.02. 

The maximum luminosity during quiescent H-burning is expressed as the sum of different terms: 

L = (18160 + 3980 log o|2)(^c -0.4468) 

_|_-|^q2.705-M.649Mc 

^;^Q0.0237(a-1.447)M^%M3=„^(l-e-^"c/o.oi) 
_-|^q3.529-(Mc,o-0.4468)AMc/0.01 

where Z denotes the metallicity, M^^ is the core mass at the first thermal pulse (see Table 2), AMc = Mc — Mcfi gives 
the actual increment of the core mass, and M^nv is the current envelope mass. Masses and luminosities are expressed 
in solar units. 

The first term (5a) represents the usual linear Mc — L relation, giving the quiescent luminosity of TP-AGB stars 
already in the full amplitude regime, with core masses in the range 0.6Mq ^ Mc ^ O.95M0. The second term (5b) 
provides a correction becoming significant for high values of the core mass, > O.95M0. The third term (5c) accounts 
for the over- luminosity produced by envelope burning (i.e. the Leb term of Eq. (2)), as a function of the envelope mass 
Mcnv A dependence on the mixing-length parameter, a, is included to reproduce the results from full calculations 
(i.e. Wagenhuber (1996) with a = 1.5, Blocker (1995) with a = 2, D'Antona & Mazzitelh (1996) with a = 2.75). 
Finally, the fourth term (5d) gives a negative correction to the luminosity in order to mimic the sub-luminous and 
steep evolution typical of the first pulses. 

The nice distinction between various terms in Eq. (5) allows us to derive the term L — Leb, that is the the 
contribution of all energy sources but for envelope burning. This is an important point indeed, since Leb is just the 
quantity we aim at evaluating by means of envelope integrations as described in Sect. 2. The luminosity produced by 
envelope burning can be eliminated from Eq. (5), by setting the corresponding term (5c) equal to unity. For the sake 
of clarity in notation, let us denote by Lm^ the resulting luminosity: 

Lm, = (18160 + 3980 log o|2)(^c -0.4468) 

_|_-|^q2.705-|-1.649Mo 
_-|^q3.529-(Mc,o-0.4468)AMc/0.01 

Hereinafter, we will refer to Eq. (6) as the standard Mc — L relation adopted in this study. 

This relation replaces the formulae by Boothroyd & Sackmann (1988a) and Ibcn & Truran (1978), for different 
ranges of the core mass, used in previous works (Marigo et al. 1996ab, 1998). The notable improvement is that the 
new prescription is based on homogeneous evolutionary calculations for a large range of core masses and metallicities. 
Moreover, it is worth remarking that the quantity L — Leb cannot be obtained when using the Iben & Truran's formula 
(1978), which is likely to mask a non-quantifiable contribution from a weak envelope burning (see Sect. 4.1 in Marigo 
et al. (1998)). 

From the above scheme it follows that the energy contribution from envelope burning, Leb, is supposed to satisfy 

the relation: 

L = Lm^ + Leb (7) 

In other words, we assume that the nuclear burning at the base of the convective envelope produces the excess of 
luminosity above the underlying Mc — L relation. 



(5a) 
(5b) 
(5c) 
(5d) 



(6a) 
(6b) 
(6c) 
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According to Wagenhuber (1996) it is possible to derive Lu from: 

log (^^^ = -0.012 - 10-l-25-113AMe _ 0.0016Menv (8) 

where the variables are expressed in solar units. Since the above relation is an analytical fit to full calculations with 
no (or quite weak) envelope burning (i.e. Xeb ~ 0), we can safely assume L = Lm^ in Eq. (8). 
Hence, we can write: 

Lh = /H^Me (9) 

where /h is the fractional contribution of radiative H-burning to Lm^ , and 

Lg + Lue - = (1 - /H)iMe (10) 

is the complementary term to Lm^, including the gravitational shrinking of the core, Lq, the shell He-burning, Lue, 
and neutrino losses, L,^. 

It turns out that, once the full amplitude regime has established, these fractional luminosities attains typical values, 
slightly dependent on the core mass, envelope mass, and metallicity, of Lu/Lm^ ~ 0.96 -=-0.98 and (Lq + Liic)/LM^ ~ 
0.04 H- 0.02. 

Finally, it is worth recalling the general validity of the method, which can be applied as well to solve the envelope 
structure of low- mass stars without envelope burning (i.e. Xeb = and L = LmJ- 

2.2. Other analytical prescriptions 

The general structure of the synthetic TP-AGB model is the same as described in Marigo et al. (1996a, 1998). With 

respect to these latter works, some input prescriptions have been updated thanks to the recent re-determinations by 
Wagenhuber (1996) and Wagenhuber & Groenewegen (1998). Besides the — L relation already quoted in Sect. 2.1, 
the other new analytical prescriptions are: 

— The core mass - interpulse period relation: 

logiip = (-3.628 0.1337 loggf^) (Mc - 1.9454) (10a) 

_-|^Q-2. 080-0. 353 log jyf2-|-0.200(Me„v-|-a-1.5) (10b) 
_-|^Q-0.626-70.30 (Mc,o-logo^) AMc ^j^q^^ 

Here three components can be distinguished, namely: the term (10a) expresses the interpulse period tip (in yr) as 
a decreasing exponential function of Mc with some dependence on the metallicity; the term (10b) gives a negative 
correction to include the effect of envelope burning in somewhat reducing the inter-pulse period; and term (10c) 
reproduces the initial increase of tip starting from values that are shorter by almost a factor of two compared to those 
derived from the asymptotic relation for the same core mass. 

— The rate of evolution of the hydrogen-exhausted core: 

^=a^ (Moyr ) (11) 
where 



q = [(1.02 ± 0.04) + O.OlTgl^] 10-" (g erg-i) (12) 
with X and Z corresponding to the hydrogen and metal abundances (mass fractions) in the envelope, respectively. 
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This work a = 2.3 BH case 




Fig. 1. Evolution of the surface luminosity as a function of the core mass for a 7Mq star with initial solar metallicity (i.e. 
Z = 0.021). The empty triangles refer to complete evolutionary calculations carried out by Blocker & Schonbcrner (1991) from 
the 1^* up to the 30**^ thermal pulse. The adopted mass-loss prescription is that suggested by Baud & Habing (1983; BH case). 
For comparison, the luminosity evolution predicted by our synthetic TP-AGB model is shown (solid line). Calculations are 
carried out starting from the same initial conditions at the first thermal pulse till the 60**" pulse. The dashed line corresponds 
to the underlying Mc — L relation adopted in this case, which is given by Eq. (6) multiplied by a factor of 1.16. The evolution 
of a 2.5Mq star with no envelope burning is also plotted. See the text for more details. 




M, (Me) 

Fig. 2. The same as in Fig. 2, but with a different prescription for mass-loss {BH4 case). The Baud & Habing's formula (1983), 
apphed since the beginning of the AGB, is then replaced with a constant rate, M = 4 x lO~'*M0yr~^, as soon the core mass 
has grown up to O.93537M0 . 
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3. Comparison with full calculations 



The most meaningful test to check the reliability of the method developed to account envelope burning is to compare 
the predictions of the synthetic TP-AGB models with those from complete AGB evolutionary calculations. We consider 
the results obtained by BS91 and B95 for a TP-AGB star with initial mass of 7Mq and chemical composition [X = 
0.739, r = 0.240, Z = 0.021]. 

We calculate the synthetic TP-AGB evolution of the 7Mq star starting from the same initial conditions at the 
first thermal pulse as in BS91, i.e. with M = 6.871M0, = O.91335M0, L = 252172.0. Specifically, the underlying 
Mc — L relation is that given by Wagenhuber & Groenewegen (1998) (Eq. (6)), multiplied by a proper factor of 1.16. 
This latter is calibrated in order to obtain the same value of the luminosity at the first thermal pulse as in BS91. 

In our envelope model the mixing-length parameter is set equal to a = 2.3 so as to obtain values of the effective 
temperature similar to those derived from full calculations (see Fig. 5 in B95), which were actually carried out with 
a lower value, a = 2.0. This is most likely due to the different opacities employed in Blocker's calculations (Cox & 
Stewart 1970), and in the present study (Iglesias & Rogers 1996; Alexander & Ferguson 1994). 

As far as mass loss by stellar winds on the AGB is concerned, two are the prescriptions here adopted as in B95, so 
that we distinguish two cases: 

1. the BH case referring to the use of the Baud & Habing's (1983) modification of Reimers' formula; 

2. the BH4 case corresponding to the use of the Baud & Habing's law (1983) from the beginning of the AGB till a 
certain stage (i.e. when Mc = O.93537M0 at the 26**^ pulse), beyond which a constant rate of 4 x 1O~*M0 yr~-^ is 
artificially introduced to mimic the onset of the superwind regime. 

For more details the reader should refer to B95. 



Table 1. A solar- met allicity 7Mq star with envelope burning: comparison between full TP-AGB modelling performed by Blocker 
& Schonberner (1991) (BS91; BH case) and our calculations (M98), for four selected values of the core mass. 



BS91/M98 


Mc 

Mq 


M 

Mq 


M 




L 






Mbol 
mag yr 




BS91 


0.91335 


6.871 


4.9 10" 


■7 


25217 




0.02 


-3.3 10" 


-5 


M98 


0.91335 


6.871 


6.8 10" 


■7 


25217 


0.000 


0.01 


-5.5 10" 


■6 


BS91 


0.92061 


6.854 


1.1 10" 


-6 


40268 




0.30 


-1.4 10~ 


-5 


M98 


0.92061 


6.850 


1.1 10" 


■6 


37788 


0.062 


0.33 


-1.5 10" 


-5 


BS91 


0.92970 


6.821 


1.7 10" 


-6 


51662 




0.46 


-7.2 10" 


■6 


M98 


0.92970 


6.807 


1.9 10" 


-6 


51537 


0.003 


0.46 


-7.6 10" 


■6 


BS91 


0.93909 


6.774 


2.1 10" 


-6 


58806 




0.52 


-4.0 10" 


■6 


M98 


0.93909 


6.745 


2.6 10" 


■6 


60760 


0.033 


0.51 


-3.4 10" 


■6 



— Mc: hydrogen-exhausted core 

— M: current stellar mass 

— M: mass-loss rate 

— L: quiescent surface luminosity 

— le^l = I ^"^lbs9i^'^^ I' Percentage difference in the evaluation of L (at given Mc) derived in the present work with respect to 
full calculations 

— fraction of the total hydrogen luminosity provided by envelope burning 

— Mboi: rate of brightening 



Figures 1 (for the BH case) and 2 (for the BH4 case) show that the results of full calculations are remarkably well 
reproduced by our synthetic calculations. The displayed luminosity evolution refers to the pre-flash maximum values 
before the occurrence of each thermal pulse. The 7M0 star starts to depart from the Mc — L as soon as it enters the 
TP-AGB phase, quickly increasing its luminosity because of the occurrence of envelope burning. 

In the BH case the luminosity is still steeply increasing when calculations have been interrupted (at the 30**^ pulse 
in BS91; at the 60*'^ pulse in this work), due to the fact that the stellar mass has not yet been significantly reduced by 
stellar winds. Differently, the artificial onset of a superwind mass-loss rate in the BH4 case causes the quick ejection 
of the envelope over the subsequent 4 — 5 interpulse periods. Our synthetic calculations quite well reproduce the 
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overall features of the luminosity evolution, comprising the initial rising, the luminosity peak, the decline following 
the activation of the high mass-loss rate, and the final re-approaching towards the Mc — L relation. 

For purpose of comparison, the predicted TP-AGB evolution of a solar-metallicity 2.5Mq star is also plotted. No 
evidence of departure from the standard Mc — L relation shows up in this case. 

As far as the 7Mq star is concerned, different properties can be then compared, namely: i) the surface luminosity 
for given core mass and envelope mass; ii) the rate of brightening Afboi; iu) the fraction of the hydrogen luminosity 
produced at the hot base of the convective envelope Leb/{Lh + Leb)', iv) the current total mass; and v) the mass- loss 
rate. These quantities are indicated in Table 1 at fixed values of the core mass for four selected thermal pulses {BH 
case; in analogy with Table 1 of BSl). 
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This work 
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BH case 
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BH4 case 



I . I I I i I I I I . I I I I . I I . i I I I I . I iH 
0.915 0.92 0.925 0.93 0.935 0.94 



Fig. 3. Percentage difference in the estimation of the quiescent surface luminosity according to synthetic AGB calculations 
(^synt) (by envelope integrations in this work (M98, solid line), and with the analytical fit of Eq. (5) by Wagenhuber & 
Groenewegen (1998)(WG98, dashed line)), with respect to the results from complete calculations (Lbs9i), as a function of the 
core mass for an evolving 7Mq star. Top and bottom panels refer to the BH case and BH4 case for mass loss, respectively. 



The agreement is indeed satisfactory. The relative errors, |ei| = |Lm98 — ^Bsgi |/-^bS91, of the surface luminosity 
estimated with our model, Lm98: with respect to that obtained in the reference work, -Lbs91: do not exceed few percents 
in all cases. The very nice accordance can be better appreciated from Fig. 3, showing (solid line) as a function of 
the core mass for the entire sequence of thermal pulses calculated by BS91 for the BH case (top panel), and B95 for 
the BH4 case (bottom panel). 

The values of cl are always quite small, also in consideration that the models are intrinsically diflPerent just because 
of the use of diflferent input physics (e.g. opacities). Moreover, it is worth noticing that in the BH case the greatest 
diflPerences in luminosity for given core mass occur just in correspondence to the first thermal pulses, the behaviour of 
which is usually irregular and model dependent, and then hardly reproducible. In the BH4 case cl is mostly confined 
to few percents as well. An increasing trend (to modest values still) shows up in the very late few intcrpulsc periods, 
which are characterised by the drastic reduction of the envelope mass. The reason of this may be ascribed to somewhat 
different core mass-interpulse period relations (in Wagenhuber's (1996) and BS91 models). Then, given the extremely 
high mass-loss rate in these stages - when envelope burning is powering down but still operating - even a small a 
difference in the core mass just prior the occurrence of a thermal pulse can correspond to rather diflferent values of the 
envelope mass, and hence of the expected surface luminosity. 
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A similar check is perforniccl using the analytical fit suggested by Wagenhuber & Groenewegen (1998; WG98) and 
quoted in Eq. (5). The term (5c), expressing the luminosity contribution from envelope burning, is evaluated by setting 
a = 2.0, as employed by BS91. The agreement with full calculations also turns out to be really good in both cases as 
illustrated in Fig. 3. 

However, the advantage of our method, based on envelope integrations during the evolutionary calculations, is that 
it not only gives an estimate of Leb and hence of the actual stellar luminosity (as the analytical fit by WG98 does), 
but the entire envelope structure of the star is consistently determined at each time step. This aspect is relevant, for 
instance, to the analysis of the nucleosynthesis due to envelope burning, which is followed in detail by integrating 
the nuclear network, once the density and temperature stratifications across the envelope are known. Moreover, the 
envelope model allows to check the sensitivity of the results to possible changes of the input physics (e.g. opacities, 
nuclear rates, mixing scheme). 
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Fig. 4. Evolution of the surface luminosity for the 7Mq AGB model as a function of the envelope mass, Menv, corresponding 
to BH4 case for mass loss. The notation is the same as in Fig. 1. Note the quick decline of the surface luminosity after the onset 
of the superwind. 



Finally, Fig. 4 shows the surface luminosity of the 7Mq star as a function of the envelope mass, with reference 
to the Bif-j case for mass loss (in analogy with Fig. 9 in B95). The agreement of synthetic results (solid line) with 
complete calculations (empty triangles) is again very good, indicating that the method is reliable to estimate the 
current strength of envelope burning as the envelope mass is being reduced by stellar winds. 

4. Synthetic evolutionary calculations 

The present paper is the first of a series devoted to present the results coming out from an extensive synthetic analysis 
of the TP-AGB evolution. 

Evolutionary calculations have been carried out to follow the TP-AGB phase (from the first thermal pulse until 
the ejection of the envelope) for a dense grid stellar models with initial masses in the range O.SM© SM©, and three 
values of the original chemical composition, [X = 0.708, Y = 0.273, Z = 0.019], [X = 0.742, Y = 0.250, Z = 0.008], and 
[X = 0.756, Y = 0.240, Z = 0.004] (Marigo 1998, PhD Thesis). For each stellar model the initial conditions at the first 
thermal pulse are extracted from full calculations performed by means of the Padua stellar evolution code (Girardi 
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& Bertelli 1998; Girardi et al. 1998, in preparation). We remind the reader that the inehision of moderate overshoot 
from core and external convection (Chiosi et al. 1992; Alongi et al. 1993) leads to a lower value of the maximum mass 
(Mup ~ SMq) for a star to pass through the AGB phase, than predicted by classical models without overshooting 

(M,p ^7-8M0). 

The analytical ingredients of the TP-AGB model are the same ones as in Marigo et al. (1996a, 1998), with some 
updates already presented in Sects. 2.1 and 2.2. We recall that the adopted prescription for mass loss is that by 
Vassiliadis & Wood (1993). The basic physical inputs employed in the static envelope model are the following. Nuclear 
reaction rates are those compiled by Caughlan & Fowler (1988), the screening factors are given by Graboske et al. 
(1973). At high temperatures (T > W^K) stellar opacities are taken from Iglesias & Rogers (1996) (OPAL), at low 
temperatures (T < 10^) the opacity tables by Alexander & Ferguson (1994) are used. The value a = 1.68 adopted in 
envelope integrations derives from the calibration of the solar model (Girardi et al. 1996), which we usually refer to as 
the standard case. For purpose of comparison, other values are also used in the computation of envelope models (i.e. 
a = 2.0 and a = 2.5). 

4.I. An improved treatment of the third dredge-up 

The usual treatment of the third dredge-up in synthetic calculations is based on the adoption of two free parameters, 
namely: the efficiency A, and the minimum core mass for convective dredge-up M™™. They are calibrated so as to fit 
some basic observational constraint, e.g. the observed luminosity function of carbon stars in the LMC (M™" = 0.58Mq 
and A ~ 0.6 — 0.7 according to Groenewegen & de Jong (1993) and Marigo et al. (1996a)). The main purpose is to 
provide useful indications on the real occurrence of convective dredge-up, given the large degree of uncertainty still 
affecting the present understanding of this process in complete analyses of thermal pulses. 

However, a weak point of synthetic models is the assumption that both dredge-up parameters are constant, re- 
gardless of the stellar mass and metallicity. On the contrary, complete models of AGB stars indicate that the onset of 
dredge-up (related to M™'") and its efficiency (related to A) is favoured in stars of higher masses and lower metallicities 
(Wood 1981; Boothroyd & Sackmann 1988b). 

In this work, this limitation of synthetic models is partially overcome. In Marigo et al. (1998, in preparation) an 
exhaustive description of the method is presented. Suffice it to formulate here the basic concepts. 

According to the results from detailed calculations of thermal pulses it turns out that, at the stage of the post- 
flash luminosity peak, the penetration of envelope convection into the inter-shell region would occur only if the base 
temperature, Tb, approaches or exceeds some critical value T^'^''^, which turns out to be almost independent from 
Mc and Z (e.g. logT^''^'^ - 6.7 as indicated by Wood (1981); logT^''"^ - 6.5 according to Boothroyd & Sackmann 
1988b). As already pointed out by Wood (1981), this useful indication provides the basic criterion to infer if and when 
dredge-up takes place. The use of M™'" is then abandoned. Of course, the critical value for T^'^'^'^ is a free parameter 
as well, which must be calibrated on the basis of some observational constraint. 

However, the improvement is real. Every time a thermal pulse is expected during calculations, envelope integrations 
are performed to check whether the condition on the base temperature is satisfied. It follows that, contrary to the 
test based on the constant Af parameter, with this scheme the response depends not only on the core mass but 
also on the current physical conditions of the envelope (i.e. the surface luminosity peak, the effective temperature, the 
mass, the chemical composition). Moreover, this method allows to determine not only the onset, but also the possible 
shut-down of dredge- up occurring when the envelope mass is significantly reduced by stellar winds. 

As already mentioned, the dredge-up parameters, A and T^'^"'^, need to be specified. The carbon star luminosity 
functions (CSLFs) in the LMC and SMC provide the observational constraint. According to Marigo (1998, Phd Thesis), 
the calibration for the LMC yields A = 0.50 and T{^""^ = 6.4, whereas for the SMC we get A 0.65 and T^""^ ^ 6.4. 
Then, a higher efficiency of convective dredge-up seems to be required at lower metallicities, which would agree with 
predictions by full calculations. 

In this way wc fix the values of the dredge-up parameters adopted during the present calculations for different 
metallicity sets. Specifically, the calibration based on the CSLF in the LMC is applied to TP-AGB models with 
metallicities Z = 0.008, the reproduction of the CSLF in the SMC specifies the dredge-up parameters for TP-AGB 
models with Z = 0.004. The former calibration is applied also to the metallicity set Z = 0.019, even if a lower A would 
be suggested by extrapolating from the Z = 0.008 and Z = 0.004 cases. 

Regarding the stellar models studied in this work (3.5M0 < M < 5.OM0), it turns out that envelope burning, 
as expected, prevents the conversion to carbon stars for most of their TP-AGB phase because of the efficient nuclear 
transmutation of newly dredged-up carbon into nitrogen. Possible transitions to the C-class may occur either early in 
the evolution when dredge-up still dominates over weak (but growing in efficiency) envelope burning, or in the very 
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late stages as envelope burning extinguishes and a few more dredge-up events possibly take place. All these aspects 
are analysed and discussed in Marigo et al. (1998, in preparation). 

4-2. Evolution in the — L diagram 

Figure 5 shows the luminosity evolution described by TP-AGB stars with initial metallicity Z = 0.008 as a fmiction 
of core mass, according to the standard case (i.e. a = 1.68). Each symbol along the curves refers to the quiescent 
pre-flash luminosity maximum, just before the occurrence of a thermal pulse. The dot-dashed line corresponds to the 
standard — L relation valid for full amplitude regime with no envelope burning. It is expressed by Eq. (6) with the 
term (6c) set equal to zero, as this refers to the first pulses which are not yet in the asymptotic regime. 

Three points are worthy to be remarked. First, all the tracks are characterised by an initial sub-luminous evolution 
typical of the first thermal pulses, below the Mc — L relation. Second, as soon as the standard relation is approached, 
stars with mass M ^ 3.5Mq increase their luminosity closely following the same relation till the end of evolution, 
regardless of the current value of the envelope mass. Quite a different behaviour characterises the luminosity evolution 
of more massive stars {3.5Mq < M < 5Mq), departing from the Mc — L relation because of the occurrence of envelope 
burning. This aspect will be discussed in detail in Sect. 5. Third, we can notice that in the cases with no envelope 
burning, stars may not exactly obey the relation shown in Fig. 5, but evolve towards slightly higher luminosities, in 
particular during the last evolutionary stages. 

This can be explained considering that the plotted Mc — L relation refers to a given value of the metallicity, 
Z = 0.008, so that the weighting factor, 3980 log(Z/0.02), in the term (6a) is constant. However, during calculations 
we take into account the possible increase of the effective metallicity, defined asZ = l — X — Y, due to the surface 
chemical enrichment produced by convective dredge-up, and possibly by envelope burning. Therefore, for a given core 
mass the luminosity is expected to be higher at increasing metallicity. The effect is more pronounced after the last 
dredge-up episodes, when the dilution of newly synthesised elements involves a smaller residual envelope mass. To 
this respect, an example is illustrated in Fig. 6, referring to the {5Mq;Z = 0.008) TP-AGB model experiencing both 
dredge-up events and envelope burning (see also Sect. 5.1). Similar results hold for the other two sets of TP-AGB 
models here considered (with metallicities Z = 0.019 and Z = 0.004). 




Fig. 5. Quiescent luminosity evolution as a function of the core mass for TP-AGB stars with initial chemical composition 
[X = 0.742, Y = 0.250, Z = 0.008], and mass as indicated nearby the corresponding track. Symbols correspond to the pre-flash 
luminosity maximum before the occurrence of each thermal pulse. The dot-dashed line represents the Mc — L relation for 
Z = 0.008. 
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5. Over-luminosity produced by envelope burning 




Fig. 6. Quiescent luminosity evolution of a 5.0Mq, Z = 0.008 star experiencing envelope burning. The open circles correspond 

to the maximum quiescent luminosity before each Ho-shcU flash. The numbers along the curve indicate the current stellar mass 
in solar units. The dashed and dotted lines represent the reference Mc — L relation for Z = 0.008 and Z = 0.016, respectively. 



As already mentioned, the energy contribution from envelope burning makes massive TP-AGB stars leave the 
Mc — L relation during part of their evolution. In this study we analyse the dependence of envelope burning on 
three basic parameters namely, the mass of the star Af, the metallicity Z of the envelope, and the mixing-length 
parameter, a = A/ Hp, vi^here A is the mixing length according to the classical theory of convection (Mixing Length 
Theory; Bohm-Vitense 1958), and Hp is the pressure scale-height. The results are clearly illustrated in Figs. 7-9, 
referring to TP-AGB stars v/ith initial mass of 3.5, 4.0, 4.5, and 5.0Mq for three values of the initial metallicity, 
Z = 0.019, Z = 0.008, and Z — 0.004, respectively. For each star of given mass and chemical composition, envelope 
integrations are carried out adopting three values of the mixing-length parameter, a = 1.68 {standard case), a = 2.00, 
and a = 2.50. 

We can notice that for particular combinations of the stellar parameters here considered (i.e. lower M, higher Z, 
and lower a, as in the (S.SM©, Z = 0.019, a = 1.68) model of Fig. 7), envelope burning does not develop or it is 
so weak that the over-luminosity effect does not show up and the the evolution of the quiescent luminosity complies 
with the Mc — L relation. Moreover, we specify that the results for the {5.0Mq,Z = 0.004) model are not actually 
present since this is just a limit case. The stellar mass slightly exceeds the critical value Mup for that metallicity, so 
that carbon ignition has occurred before the onset of the thermally pulsing regime. 

It is worth remarking that M, Z, and a significantly affect the luminosity evolution of a star with envelope burning 
as long as the stellar mass is not so dramatically lowered, before the re- approaching onto the Mc — L relation. This 
corresponds to most of the TP-AGB duration, since high mass-loss rates are typically attained at the very end of the 
evolution (Vassiliadis & Wood 1993). Once the super- wind regime has developed, the subsequent luminosity evolution 
is crucially controlled by the efficiency of mass loss in reducing the envelope mass down to the complete ejection. 
An earlier onset of the super-wind would favour an earlier extinction of envelope burning and, correspondingly, a 
lower maximum luminosity would be reached by the star before the subsequent decline towards the Mc — L relation. 
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In principle, the onset of the super-wind regime as soon as the star enters the AGB phase may even prevent the 
occurrence of envelope burning. 

In the following discussion we will consider the effects produced by varying M, Z, and a, for a given mass-loss 
prescription (i.e. Vassiliadis & Wood 1993). The main results relevant to the present analysis are indicated in Table 2. 

5.1. The dependence on M 

For fixed Z and a, envelope burning is expected to be more efficient in stars of greater mass. From Figs. 7 - 9 it 

turns out that the maximum excursion of the luminosity above the Mc — L relation is higher in initially more massive 
models (see also the M^^f^ enter in Table 2). Moreover, the development of envelope burning crucially depends 
on the current envelope mass during the evolution. This effect is exemplified in Fig. 6, referring to the case of the 
{5Mq,Z = 0.008, a = 1.68) model. This behaviour of the luminosity agrees with the results from full evolutionary 
calculations of the TP-AGB phase (e.g. B95; Boothroyd & Sackmann 1992; Vassiliadis & Wood 1993). 

Envelope burning develops since the first sub-luminous interpulse periods, so that when thermal pulses attain the 
full-amplitude regime the star does not settle on the Mc — L relation (dashed line), but quickly reaches higher and 
higher luminosities, with a brightening rate (few times —10"^-; — 10~^ mag yr~^; see the Mboi entry in Table 2), much 
greater than expected from the slope of the standard relation (typically —7-. — 9 x 10~^ mag yr~^). 

The luminosity growth goes on as long as the envelope remains massive enough to support high envelope base 
temperatures (Tb > 40 — 60 x 10^ K). In fact, the luminosity decline after the maximum is concomitant with the 
onset of the super-wind, when the star starts to rapidly lose mass at significant rates (i.e. M ~ 10~^ 1O~*M0 yr~^). 
This drastically weakens the efficiency of envelope burning so that, finally, the over- luminosity vanishes and the star 
approaches again the — L relation where it remains during the very last stages till the end of evolution. To this 
respect, we can notice that the final recovering of the — L relation (dotted line) is consistent with an effective 
metallicity {Z = 0.016) that is twice the value at the beginning of the TP-AGB phase {Z = 0.008; dashed line), 
because of the occurrence of both dredge-up and envelope burning. 

In Table 2 we also indicate the value of the envelope mass at the shut-down of envelope burning (Menv(noEB) 
entry), i.e. when its contribution to the stellar luminosity has decreased to less than 1 %. Note that Monv(noEB) is not 
constant at all, but it varies mostly in the range between roughly 2Mq and 0.5Mq. It results that the more efficient 
envelope burning has been (e.g. for larger M, lower Z, higher a), the smaller is this critical value of the envelope mass. 

5.2. The dependence on Z 

At given M and a, the over-luminosity is more pronounced at lower metallicities. It is worth remarking that the direct 
effect of chemical composition on the base temperature is negligible itself (Sackmann & Boothroyd 1991). The greater 
efficiency of envelope burning is mostly due to the fact that a star of given M enters the AGB phase with a core 
mass, Mc, that is greater for decreasing metallicity (see for example the Mc,o entry in Table 2). Considering that in 
the deepest layers of the envelope close to the core, the radiative gradient depends like Vr oc M^, where 7 is always 
positive (0.5 2.5, as indicated by Scalo et al. 1975), it follows that a deeper penetration of external convection is 
expected for higher values of the core mass, and hence for lower metallicities. 

5.3. The dependence on a 

The dependence of envelope burning on a is remarkable as already pointed out by various authors (Sackmann & 
Boothroyd 1991; D'Antona & Mazzitelli 1996). An increase of a causes an overall local rise of the temperature profile 
across the envelope, so that both the effective temperature, Tcff, and the base temperature, Tb, are hotter. It follows 
that a affects both the external configuration of an AGB star (i.e. the position on the H-R diagram), and its internal 
structure (i.e. the temperature stratification of the convective envelope and related nuclear energy generation). In 
particular, the higher temperatures, attained at the base of the convective envelope at increasing a, determine a 
stronger efficiency of nuclear burning. This is evident in Figs. 7-9 from the greater amplitude of the luminosity 
excursion in models with the same mass and metallicity (see also the M^'^ and (iEB/-t')max entries in Table 2). 

However, changing the value of the mixing-length parameter produces important consequences which are not merely 
related to the efficiency of envelope burning, but deal with further aspects of the evolution of these stars. 

A notable point concerns mass loss, and consequently the TP-AGB lifetimes (see the ttp-agb entry in Table 2). 
From our analysis it turns out that for a given model, the onset of the super-wind regime is delayed for higher values 
of the mixing-length parameter. Consequently, for given initial mass and metallicity the whole duration of the TP- 
AGB phase is longer for a star with a more efficient envelope burning obtained by increasing a. This result does not 
contradict the result that the high luminosities produced by envelope burning would anticipate the end the evolution 
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Table 2. Properties of TP-AGB stars with envelope burning. 
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— Z: initial motallicity 

— Mi: initial mass at the ZAMS [Mq) 

— Mcfi- core mass at the onset of the TP-AGB phase (Mq) 

— -Ttp-agb: TP-AGB lifetime (yr) 

— Mf: final mass (Mq) 

— M^'j^: bolomotric magnitude at the tip of the AGB 

— M^^j^'^: bolometric magnitude at the maximum efficiency of envelope burning, before the onset of the super- wind phase 

— (^EB/ijinax: maximum relative contribution of envelope burning to the surface luminosity 

— Menv(noEB): envelope mass at the extinction of envelope burning (Mq) 

by triggering enhanced mass-loss rates (see for example BS91). To clarify this point, let us consider the competition 
between luminosity and effective temperature in determining the efhciency of mass loss. According to Vassiliadis & 
Wood's prescription (1993), adopted in our calculations, the dependence of the mass-loss rate before the development 
of the superwind, can be expressed as log M oc P oc T^'^^ L'^-'^'^ M^^-^ , where P is the fundamental period of pulsation. 
Hence, an increase of a operates in two opposite directions. From one side, it favours mass loss via the term L^'^'^ , as it 
strengthens the efficiency of envelope burning, quickly leading to higher luminosities. From the other side, it weakens 
mass loss via the term T^^'^^ , as it reduces the stellar radius, yielding higher values of the effective temperature. The 
net result depends on the prevailing effect, that in our case turns out to be related to the increase of the effective 
temperature. It is interesting to notice, however, that different formulations for mass loss could produce different 
results, owing to their specific dependence on stellar parameters, i.e. luminosity and effective temperature. 



14 P.Marigo: Envelope burning over-luminosity: a challenge to synthetic TP-AGB models 



80 



60 



© 

CO 

O 



40 



20 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

-M = 3.5 Mq 

;z=o.oi9 


1 1 1 1 1 1 1 1 


- — a=1.68 
- - a = 2.00 
■ a=2.50 


/ 
/ - 
/ - 

/ 


" .-^^ 
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 



80 



60 - 



© 

O 



40 - 



20 



I I I I I I I I I I I I I I I I I I I I I I I I I I 
"M=4.0 



Z = 0.019 



0.7 0.8 0.9 1 1.1 1.2 
Me (Me) 



I I I I I I I I I I I I I I I I I I I I I I I I I 
0.7 0.8 0.9 1 1.1 1.2 



Me (Me) 



80 



60 - 



© 

CO 

O 



40 - 



20 



I I I I I I I I I I I I I I I I I I I I I I I I I 
■M=4.5 Mq 
"Z = 0.019 



I I I I I I I I I I I I I I I I I I I I I I I I I I 



80 



60 - 



© 

CO 

O 



40 - 



20 - 



I I I I I I I I I I I I I I I I I I I I I I I I I I 
"M = 5.0 



Z=0.019 




I I I I I I I I I I I I I I I I I I I I I I I I I I 



0.7 0.8 0.9 1 1.1 1.2 
Me (MJ 



0.7 0.8 0.9 1 1.1 1.2 

Me (MJ 



Fig. 7. Luminosity evolution of 3.5, 4.0, 4.5, 5.0 Mq TP-AGB model stars with initial metallicity Z = 0.019 as a function of the 
core mass for three values of the mixing-length parameter, a, as indicated. Stellar luminosities refer to the pre-flash maocimum 
stage. The long-dashed lines correspond to the reference standard core mass-luminosity relation. 
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Fig. 8. The same as in Fig. 7, but with Z = 0.008. 
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Fig. 9. The same as in Fig. 7, but with Z = 0.004. 
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6. Concluding remarks 

In this study it has been shown that the break-down of the — L relation and the related over-luminosity effect 
produced by the occurrence of envelope burning in the most massive TP-AGB stars can be consistently taken into 
account in synthetic calculations. In fact, the original method based on envelope integrations is able to successfully 
reproduce the results of complete evolutionary models. 

This nice result strengthens the reliability of envelope integrations in synthetic TP-AGB models for various related 
aspects, such as the nucleosynthesis occurring in hot bottom envelopes and the chemical yields. 

The results of the present paper derive from extensive calculations of the AGB phase carried out for a fine grid of 
stellar masses (O.SMg <, M < 5Mq) and various metaUicities {Z = 0.019, Z = 0.008, and Z = 0.004). This work is 
actually a part of a project aimed at investigating various aspects of the AGB evolution by means of a flexible and 
accurate synthetic model. If the flexibility is an intrinsic characteristic of synthetic codes, this model gains in accuracy 
thanks to the use of analytical prescriptions derived from detailed AGB calculations (e.g. Wagcnhuber 1996), coupled 
to a complete envelope model updated with recent input physics (i.e. opacities, nuclear reactions rates). 

Up-coming papers are devoted to present and analyse other issues of interest. In particular, basing on the improved 
treatment of the third dredge-up included in synthetic calculations (briefly outlined in Sect. 4.1) we will address the 
question of reproducing the luminosity functions of carbon stars in both the LMC and SMC (Marigo et al. 1998, in 
preparation). A further work will be dedicated to present the results on the predicted changes in the surface chemical 
abundances of AGB stars and related chemical yields, as a function of the stellar mass and metallicity (Marigo et al. 
1998, in preparation). 
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